Anisotropy on the Fermi Surface of the Two-Dimensional Hubbard Model 
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We investigate anisotropic charge fluctuations in the two-dimensional Hubbard model at half 
filling. By the quantum Monte Carlo method, we calculate a momentum-resolved charge com- 
pressibility K(fc) = d{n{k)) /dfi, which shows effects of an infinitesimal doping. At the temperature 
T ~ t^/U, K{k) shows peak structure at the (±7r/2, ±7r/2) points along the \kx\ -\- \ky\ = tt line. A 
similar peak structure is reproduced in the mean-filed calculation for the d-wave pairing state or the 
staggered flux state. 
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Effects of electron-electron interaction have attracted 
lots of interest. One of the most significant phenom- 
ena in the strongly correlated electron system is the 
Mott transition, which is a quantum phase transition 
driven by the interaction. This Mott insulator has a 
finite charge gap and the anitiferromagnetism often ac- 
companies. Recently, various anomalous properties have 
been found near the Mott transition. Especially, there 
has been a proposal that the interaction brings about 
anisotropy in the low energy excitations. For exam- 
ple, the spectral weight of the-twQ-dimensional Hubbard 
model has been investigatedBo^Ba. In addition, the de- 
formation of the Fermi sucface due to the interaction 
occurs in the t — J modelHQ. The singular momentum 
dependence is also observed experimentally. Angle re- 
solved photoemission spectroscopy (ARPES) measure- 
ments suggest anisotropic properties in the low energy 
excitationsQ. These anomalous properties are widely ob- 
served in the strongly correlated electron system and they 
could be an evidence of the non-Fermi liquid behavior. 
It is thus an important and appealing topic to study the 
low energy excitations with momentum resolution. 

The anisotropy in the low energy excitations is nat- 
ural, rif we assume the c?-wave pairing state or the flux 
statda. They both have a singular energy dispersion that 
has gap nodes along the diagonal directions \k,j.\ = \ky\, 
while a gap opens around the (±7r, 0) and (0, ±7r) points. 
Thus, one expects anisotropic charge excitations near the 
Fermi surface. In general, however, these states compete 
with other instabilities. Especially, at a rational filling, 
the system often belongs to the antiferromagnetic Mott 
insulator, where the charge degree of freedom is frozen. 
In that case, the Neel state well describes the ground 
state. Nevertheless, the d-wave pairing state or the flux 
state can give a sound basis for the interpretation of some 
singular phenomena near the Mott insulator at a finite 
temperature. 

In this letter, we investigate the charge fluctuation 
in the half-filled Hubbard model on a two-dimensional 
square lattice. Using the quantum Monte Carlo method, 
we calculate the momentum-resolved charge compress- 



ibility 



d{n{k)) 



where {n{k)) = (c^Cfe) denotes the momentum distri- 
bution function and fi the chemical potential. The in- 
tegral of K{k) is equal to the charge compressibility 
K = J dkK{k). If system has a finite charge gap, the 
charge compressibility decreases exponentially towards 
zero as a temperature is lowered. Thus, at half filling, 
we cannot calculate K(fc) with sufficient numerical accu- 
racy at a very low temperature. At a finite tempera- 
ture, an infinitesimal shift of /i causes an infinitesimal 
doping. This means K[k) reflects the distribution of in- 
finitesimally doped carriers with momentum resolution 
5n{k) « K{k)Sfj,. Without the interaction, K{k) is peaked 
on the Fermi surface and its value is constant. At half 
filling, this non-interacting Fermi surface is a square in 
the Brillouin zone, \kx \ + \ky\ = tt. We define the Fermi 
surface(FS) as this square in this paper. 

We calculate K{k) at half filling. Here the charge 
degree of freedom is almost frozen and the system is 
dominated only by the insulating fixed point (antiferro- 
magnetic Mott insulator) at the low temperature (T ^ 
t^/U). On the other hand, at a sufhciently high 
temperature(T ^ J7), the Coulomb interaction U is irrel- 
evant. We focus on the intermediate temperature region 
(T /U), expecting an interaction between charge and 
spin degree of freedom gives non-trivial feature on K(k) 
even at half-filling. 

The Hamiltonian of the two-dimensional Hubbard 
model is given by 

+ (n,| - 1/2) (n,i - 1/2) - m^".-, 

i i.a 

where denotes the nearest-neighbor links, and t the 
nearest-neighbor hopping amplitude. The system is on 
a square lattice and we impose periodic boundary condi- 
tions. 
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In order to obtain approximation-free results, we 
employ the finite temperaturfe. . auxi liary field quantum 
Monte Carlo (QMC) methodEl'By. In this method, 
physical observables are evaluated in the grand canon- 
ical ensemble. This makes it possible to obtain K(fc) by 
direct sampling in the QMC simulations as 

K{k) = /3 ( Hfc)^n(fc)) - (n(fc))(^n(fc)) ) , 

V fe k / 

where /3 denotes an inverse temperature. If we calcu- 
late K{k) by numerically differentiating the QMC data 
(n(fc)), we need to perform simulations for the doped 
system. This brings about the notorious sign problem, 
which prevents us from getting reliable data. On the 
other hand, the direct evaluation of K{k) only needs a 
simulation at half filling, where the sign problem does 
not occur due to the particle-hole symmetry. The num- 
ber of the electrons is not fixed at half filling in the QMC 
ensembles, even if we set the chemical potential fi — 0. 
Thus, the information on infinitesimally doped systems 
is statistically taken into account in /t(fc). 

The simulations were performed on a 16 x 16 square 
lattice. The finite size effects on K{k) are not observed 
with this lattice size at the temperatures we have studied. 
The Trotter time slice size is set to be At ~ 0.1/t. We 
have checked that the systematic error due to the Trotter 
decomposition does not change qualitative features. For 
the interacting case, the strength of the interaction is set 
to be U/t = 4, where the charge gap is Eg ~ 0.6ni3. 
We have typically performed 500 Monte Carlo sweeps in 
order to reach a thermal equilibrium followed by 10"* 
measurement sweeps. The measurements are divided into 
10 blocks and the statistical error is estimated by the 
variance among the blocks. 

At first, let us discuss results at the temperature 
T - t^/U. The results of K{k) for U/t = and U/t = 4 
are shown in Fig. ^. For the non-interacting case, the 
value of K{k) on the FS is constant. On the other hand, 
for U/t = 4, peak structure emerges in K(fe) at the 
(±7r/2,±7r/2) points. It indicates that the (±7r/2,±7r/2) 
points are more sensitive to the shift of the chemical po- 
tential than (±7r,0) or (0,±7r). In other words, the sys- 
tem is more compressible at the (±7r/2, ±7r/2) points. 
Here we note that the interaction does not change the 
line defined by {n{k)) = 0.5, which is identical with the 
FS in the non-interacting case. In that sense, the shgype 
of the Fermi surface itself is not deformed at half fiUingliS. 

Next, we discuss the temperature dependence of K{k). 
The results of K{k) for U/t = 4 at several tempera- 
tures are provided in Fig. g. The peak structure at the 
(±7r/2, ±7r/2) points is clearly observed at T ~ t^ /U . 
It becomes ambiguous as the temperature increases and 
vanishes above T ^ U . Here we note that the antifer- 
romagnetic correlation length is smaller than the linear 
system size at these temperatures {T/t > 0.2). There are 
two known characteristic energy scales in the Hubbard 
model. One is the Coulomb interaction U and the other 



is the effective superexchange interaction J ~ t'^/U. Our 
results show the peak structure in K{k) emerges at the 
temperature T ^ J and vanishes at T > [/. 

In the half-filled Hubbard model, an antiferromagnetic 
long-range order appears at T = 0. At a finite tem- 
perature, the system does not have any long-range or- 
der. Then there is no a priori reason to expect such 
an anisotropy in the charge compressibility. Thus, let 
us compare the QMC results with various mean-field so- 
lutions at the temperature T ~ J. We focus on three 
possible mean-field solutions; the Neel state, the d-wave 
pairing state, and the staggered flux state. As an effec- 
tive model of the Hubbard modelliJ, we use the following 
t — J model to compare with our QMC results 

{j.k).a 




where Si — | X^o- o-' •^L'^'o-o-'Cjcr' and the double occu- 
pancy at the same site is prohibited. We take order pa- 
rameters as 

^jk = -^((cj-TCfci - CjiCfct)) (d-wave pairing), 

X]k = ^((cjTcIt + '^ii'^li) ) (staggered fiux), 
rrij = J(5|) (Neel). 

We set Ajfc to have the d-wave symmetry and its am- 
plitude is constant A. The link order Xjk is chosen so 
that the effective magnetic fiux for each plaquette takes 
(f) or —0 alternately. The parameters satisfy A/i < 1 or 
(f) < n where the low-lying-excitations are described by 
anisotropic Dirac fermionsE£l. The staggered magnetiza- 
tion TTij is (— l)^^+-'"m. In order to apply the mean-field 
ansatz, we introduce auxiliary fields to decouple the su- 
perexchange term and incorporate the constraints of no 
double occupancy. Next, we apply the saddle-point ap- 
proximation and neglect the fluctuation around it by tak- 
ing order parameters as an input parameter (without self 
consistency condition). In low dimensions, the fluctua- 
tions generally play a role of destructing the long-range 
order and it is natural to expect it in our case. Our focus 
is, however, not on such a long-distance behavior. 

We show K{k) at half filling for these mean-field states 
in Fig. ^. The peak structure at the (±Tr/2,±7r/2) points 
is observed in the d-wave pairing or the staggered flux 
state, which is similar to the QMC results. This implies 
the d-wave pairing state or the staggered flux state com- 
petes with the Neel order at a finite temperature and 
they give a good description of the short-distance be- 
havior. However, we don't claim that these mean-filed 
states become long-range ordered at half filling. Our in- 
terest is in which type of the order parameter can give 
the anisotropy in the charge compressibility at T ~ J. 
Even if these mean-filed states are unstable, they can 
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still play an important role at a finite temperature. On 
the other hand, K{k) of the Neel state gives a constant 
value on the FS, which does not reproduce the QMC re- 
sults. The ground state of the half-filled system is well 
described by the Neel state, while the system behaves 
essentially as a non-interacting case at a sufficiently high 
temperature {T U). Both the Neel state and the non- 
interacting metallic state give a constant K(fc) on the FS. 
Therefore, one of the natural scenarios for the anisotropy 
in nik) is that different kinds of fixed points exist and 
bring about some singular phenomena at an intermediate 
temperature. Finally, we note that a possible admixture 
of different orders may happen at a low temperature. 

In summary, we have investigated the charge ffuc- 
tuation in the two-dimensional Hubbard model by the 
quantum Monte Carlo method. The momentum-resolved 
charge compressibility K(k) is focused on at a finite tem- 
perature. It gives information on the infinitesimal dop- 
ing in the Mott insulator. The peak structure at the 
(±7r/2,±7r/2) points is observed in K{k) at the temper- 



ature T ^ /U . It is qualitatively consistent with the 
calculation for the d-wave pairing state or the staggered 
fiux state, while the low temperature(T ^ /U) behav- 
ior is dominated by the antiferromagnetic Mott insulator. 
This peak structure disappears at the high temperature 
T > U where the Coulomb interaction is irrelevant. The 
crossover observed in our results refiects the existence 
of several fixed points (including unstable ones) in the 
strongly correlated electron systems. 
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FIG. 1: The momentum-resolved compressibility K,{k) — 
d{n{k)) / dfi obtained by the quantum Monte Carlo simula- 
tions at T/t = 0.2 on a 16 X 16 lattice: (a) for U/t = 0; (b) for 
U/t = 4. The areas with large values of K(fc) are highlighted 
for the contour plots. Without the interaction((7/t = 0), 
K(fe) is constant on the Fermi surface. On the other hand, 
for U/t — 4, K(fe) shows peak structure at the (±7r/2, ±7r/2) 
points. 
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FIG. 2: The momentum-resolved compressibility K(fc) = 
d{n{k)) / obtained by the quantum Monte Carlo simula- 
tions on a 16 X 16 lattice for U /t = 4 a.t various temperatures. 
The peak structure at the (±7r/2, ±7r/2) points vanishes at 
T>U. 
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FIG. 3: The mean-field calculations of the momentum- 
resolved compressibility K{k) = d{n{k)) / djj. at T/t = 0.2: 
(a) for the d-wave pairing state(A/t = 0.2); (b) for the stag- 
gered flux state {(p = 7r/6); (c) for the Neel state(m/t = 0.6). 
The areas with large values of K{k) arc highlighted. The K{k) 
shows peak structure at the (±7r/2, ±7r/2) points on the Fermi 
surface for the d-wave or staggered flux state, which is con- 
sistent with the quantum Monte Carlo simulations. On the 
other hand, for the Neel state, K{k) is constant on the Fermi 
surface. 



